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1 INTRODUCTION 


This paper describes a Fortran IV simulation study of the all-digital phase- 
processing circuitry designed by the NASA Omega staff of the Ohio University 
Avionics Engineering Center, A digital phase-lock loop (DPLL) forms the heart 
of the Omega navigation receiver prototype. Through the DPLL/ the phase of the 
10.2 KHz Omega signal is estimated when the true signal phase is contaminated 
with noise. The study has provided a convenient means of evaluating loop per- 
formance in a variety of noise environments/ and has proved to be a useful tool for 
evaluating design changes. 

The DPLL uses a frequency synthesizer as the reference oscillator. 

The synthesizer is composed of synchronous rate multipliers (SRM’s) driven by a 
temperature-compensated crystal oscillator (TCXO). Use of the SRM's in this 
application introduces phase {itter which degrades system performance. Simulation 
of the frequency synthesizer is discussed in Part III of this. paper. 

The goals of the simulation were fourfold: 

A, To analyze the circuits on a bit-by-bit level in order to evaluate the 
overall design; 

B, To see easily the effects of proposed design changes prior to actual 
breadboarding; 

C. To determine the optimum integration time for the DPLL in an environ- 
ment typical of general aviation conditions; and 

D. To quantify the phase error introduced by the SRM synthesizer and' 
examine its effect on the system. 



II THE DIGITAL PHASE-LOCK LOOP 


A. Introduction. DPLL's have been the subject of intense study in the past 
few years (^/^/^). Application of phase-lock techniques to digital systems has 

led to loop filters of first, second, and higher orders. The DPLL described in this 
section is a first-order loop which stores information about the relative phase dif- 
ference between the local oscillator and the Input signal . 

B. Description of Circuit. The DPLL follows closely the ideas of 
J.M. Clark^^^, first presented in 1968. The DPLL presently in use at Ohio 
University is essentially a highly adapted version of his digital phase tracking 
filter. A detailed discussion of the loop circuit operation is given by Chamberlin^^^ . 

The circuit employs a 6-bit counter which cycles once every 64 pulses 
of the 652.8 KHz clock (see Figure 1). This counter therefore cycles at a rate of 
10.2 KHz (=652.8 KHz f 64). The number in the counter is compared continuously 
to the six most significant bits in the up-down counter, which is in turn fed by the 
phase detector. During the brief period when the numbers in the two counters are 
bit-wise identical, the output of the comparator goes to the logical 1 state. Al- 
though the duration of this pulse may vary, the uni -directional counter generally 
advances at a much faster rate, indicating a pulse width of 1/652.8 KHz, or 
approximately 1 .5 ^sec. This pulse is used to trigger a monostable, creating a 
window whenever the contents of the two registers are identical . 

When viewed through this window, the zero crossings of the Incoming 
signal are seen to be either leading or lagging in relative phase. When a lag is 
indicated, the up-down counter receives a down-count command, and the next 
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window is created earlier in time, in this way the edge of the window gating 
function is brought more and more closely into coincidence with the zero crossings 
of the true Omega signal. Similarly, when a phase lead is indicated, count-up 
commands cause the window to be created later in time. 

When lock is obtained, the six most-significant bits in the bi-directional 
counter represent the relative phase difference between the input frequency and 
the local oscillator (LO). The LO can be considered as a 10.2 KHz clock, whose 
zero crossings coincide in time with the zero-count of the 6-bit counter. 

C. The DPLL as a Cross-Correlation Device. In a cross-correlation type 
receiver scheme, two signals are combined to give a d.c. level by first delaying 
one signal in time, multiplying their resulting signals, and integrating or iowpass 
filtering (see Figure 2). To see how the DPLL performs this same function in an 
all-digital sense, it is first necessary to view the circuit operation under no-noise 
conditions. Under these conditions there is no phase shift due to noise perturbations, 
nor is there a Doppler shift. The only phase error is then due to initial phase dif- 
ference at turn-on . 

Assume that the initial phase difference is such that a count-up command 

is given to the up-down counter. If N bits are being integrated, there will be a 

delay of 2^ up-count commands before the comparator logic "sees" a difference 

in time between successive occurrences of identical numbers in the tv/o counters. 

N 3 

The time delay in such a case is (2 /10.2 x 10 ) seconds. This integration 
performs the low-pass filtering operation of h(t) in Figure 2. 

The output of the comparator, ©, is the estimate of the Omega signal 
phase, When lock is obtained, and the binary number contained 
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Figure]. The DPLL Circuit. 



Figure 2. A Generalized Cross-Correlation Device. 
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in the six most significant bits of the up-down counter is the phase difference, 
analogous to the d.c.ievel obtained in an analog cross-correlation device, 

D. Description of the Simulation, The digital nature of the circuit 
indicated the use of the digital computer as the most convenient simulation tool . 
Since the circuit was already built and was, in fact, yielding phase information' 
which compared favorably with the Tracor series 599R receiver used for control- 
group data collection purposes, the circuit model was wholly deterministic. 

The generation of zero-crossing information for the Omega signal in 
noise, however, was necessarily stochastic in nature, A complete description of 
the signal -generating procedures is given in Part I! E. 

CORDET is based on the technique of periodic scanning; that is, each 
portion of the DPLL is observed and updated with each clock pulse, significant 
information about the state of the system is recorded, and the process is repeated. 
In the case of CORDET, the significant information includes the states of all logic 
circuit elements, including phase difference information referenced to a local 
oscillator. The phase information is, of course, of primary importance in an 
Omega navigation receiver. The ability to compare the DPLL output with the true 
phase of the simulated Omega signal, and to provide graphical output, are ad- 
vantages of a computer simulation over a laboratory hardware experiment. A 
user’s technical description is provided in Appendix A. 

E. Signal input Simulation, Since the receiver front. end provides a 
15-Hz bandwidth narrow-band filter, the derived phase information can be con- 
sidered to result from a noise-plus-signal input to a narrow-band system. 

The composite input signal can be expressed as 
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0 ) 


g(t) = f(t) + A cos {c^^i + r ) 

where is the cenfer frequency of the filter, and f(t) is the narrow-band 
noise process 

f(t) = E(t) cos (w^t + .}!^(t)). (2) 

(3) 

It can be shown^ ' that the amplitude process E(t) is Rayleigh distributed, 
while the phase process j^(t) is uniformly distrubuted on the interval (0,2-n'). 

Now, f(t) can be written in the form defining its quadrature com- 
ponents, 

f(t) = N^(t) coso)^t + N2(t) sinWj^t (3) 

(3) 

where N](t) and N 2 (t) are independent Gaussian processes' , 

From equation 1, g(t) may be expressed as 

g(t) = f(t) A cosy’costo^t + A sin^^sinw^t . (4) 

Substituting for f(t), and using 

cos (a + b) = cos a cos b - sin a sin b, 

we have 


g(t) = (N^ (t) + A cos'fr) cosw^t + (N 2 (t) + A sin^) sincOj^t. 
Equation 5 may be put in the form 

g(t) = R(t) cos(w^t + ©(f)) 

where 

R(t) = |^[N'](0 A cos^ + jN2(t) + A sinjjrJ^^ ^ 

and 


N 2 (t) + A sin jir 

©(t) = arctan 

N^ (t) + A cos ~jr 


(5) 

( 6 ) 

(7) 

( 8 ) 
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Equation 8 is therefore the input phase representation to the phase 
detector, where (t) and N 2 (t) are statistically-independent Gaussian white 
noise processes with zero mean and unit variance, and A is the normalized signal 
amplitude, in general can be a function of time to represent the Doppler 
shift in the dynamic case. At speeds less than Mach 1, however, the Doppler 
shift is well below the 5,06° (1/64 of a 10.2 KHz cycle) quantizing interval of • 
the DPLL, and so be considered constant over the sampling interval . 

A block diagram of the noise input phase simulation is given in 

Figure 3. 

F. Conclusions. In determining the optima! integration time of the first 
order loop, consider the case of a small aircraft in a noisy environment. As 
the time-multiplexed Omega transmission sequence progresses, the aircraft will 
change position between successive transmitting time Intervals of the 10.2 KHz 
frequency. For an aircraft flying a course coincident with a station pair base- 
line (worst case), the change in position will amount io 12° of 10.2 KHz phase 
at a speed of 200 mph, conservative for private aircraft. For this reason, step 
phase inputs of 12° were taken both in a noise environment of 10 decibels SNR 
and in the no -noise case. 

To determine ojatimal integration time, phase error at the end of the 
625 msec, sampling interval is plotted vs. integration time in Figures 4 and 5. 

As could be expected, the best integration time in noise is one which is as long 
as possible without leaving a residual phase error at the end of a sampling Interval 
in the no-noise case. For aircraft at 200 mph, 0.1 sec., or 10 bits of integration 
proves to be optimal . 


- 7 - 



costot 



N cosut 

2 
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Figure 3, Block Diagram Simulation of Narrow-Band Noise. 
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Figure 5. 


Bits of Integration 

Typical Plot of Phase Error at End of Sampling Interval 
vs. Bits of Integration, 10 Decibel SNR. 
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Ill SRM FREQUENCY SYNTHESIZER 


A. Introduction, Many applications in communications and navigation 
require synchronous signal detection using a local oscillator with high stability, 
and at non-standard frequencies. For example, most Omega navigation receivers 
operating in hyperbolic mode require that a local oscillator (LO) at some integral 
multiple frequency of 40.8 KHz be used for detection. The 40.8 KHz frequency 
is the least common multiple of 10.2 KHz and 13.6 KHz (I.e., 4 x 10.2 and 
3 X 13.6). These are two of the three frequencies broadcast by the Omega system, 
the other being 1 1 .333 KHz. 

The Omega receiver prototype under development at Ohio University 

^81 

uses a LO frequency of 2.6112 MHz (64 x 40.8 KHz) . It is required that 

7 

the LO have a long-term stability of at least one part in 10 . To meet this 
requirement the LO must be of the temperature-compensated crystal oscillator 
type (TCXO), Since 2.6112 MHz is not a standard -frequency crystal, and 
special cutting of the crystal would be expensive, if was decided to synthesize 
the required frequency from a 5 MHz TCXO standard. The synthesized 2.6112 MHz 
frequency is then divided by 4 and the resulting 652.8 KHz wave (64 x 10.2 KHz) 
is used as the reference in the digital phase-lock loop described in Part II A. 

Use of the 5 MHz TCXO reduces the problem of long-term drift in the 
LO, but a new short-term phase jitter is introduced by the frequency synthesizer. 
This phase fitter, seen as additive noise when referred to the input signal off 
the air, must first be quantified before system degradation and design improvements 
can be evaluated. 
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8. The SRM Synthesizer. 


Rate multipliers have been in use in discrete 


(9) 

form for several years ^ and recently (1972) programmable versions have become 
available in integrated -circuit form^^^^. Rate multipliers are essentially counters 
which count moduio-N. When supplied with appropriate logic circuitry, the 
counter will gate clock pulses in a predetermined sequence, effectively, multiplying 
the input clock rate by a factor of M/N, v/here M is determined by the programming 
circuitry. 

Literature concerning synchronous rate multipliers (SRM's) appears to be 
scarce, although the design employed by the Avionics Research Center staff at 
Ohio was inspired by references 9,10, and 11 . The synthesizer utilizes the 
SN 741 67 Decimal Rate Multiplier (DRM) manufactured by Texas Instruments, 

Each DRM multiplies the input frequency by M/lO, where M can be programmed 
to any integer from 0 through 10. 

A block diagram of the frequency synthesizer is shown in Figure 6. 

It consists of a 5 MHz TCXO followed by a chain of DRM's, The DRM chain 
is programmed in binary to multiply the input-rate by .52224. That is, when- 
ever 10 clock pulses enter the chain, 52,224 pulses are produced at the 
OR-gate output. To be able to calculate phase jitter in the output, it is clear 
that the sequence of pulses produced by each DRM must first be known. A 
FORTRAN computer program was written to determine the DRM outputs. 

C, DRM Simulation. Each DRM is composed of four T~type flip-flops 
and accompanying control logic gates (see Figure 7), Using the truth table 
associated with each of the logic elements, a program was written for the IBM 
System 360 computer using FORTRAN logical variables. The program, which 
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BINARY-CODED DECIMAL RATE MULTIPLIERS 
5 2 2 2 4 



FIguce 6, DRM Frequenc/ Synthesizer. 










Ennhifa 



Figure 7. Funcfional Block Diagram of SN 54167/SN 74167 Synchronous Decade Rate Multiplier. 





takes 10 seconds to execute, is listed in Appendix B. 

The results of the computer mode! are summarized in Figure 8» From the 
timing diagrams it is evident that programming the DRM inputs serves only to 
delete certain of the clock pulses. When the DRM's are cascaded in order to 
multiply by a longer decimal fraction, pulses are added whenever the ENABLE 
line to lower-order stages goes LOW. As seen from the ENABLE output in 
Figure 8, this con occur only for a period of one clock pulse in every ten. 

Note that the outputs determined by the simulation are identical to those 
produced by the circuit in the laboratory. Programmings of 5, 2, and 4 are shown 
in Figure 9. Results for all programmed inputs of Figure 8 were in agreement 
with laboratory observations. 

it is also evident that the first stage in the DRM chain will determine the 
worst case for phase jitter. For example, when the DRM is programmed to mul- 
tiply by .5, the output will be the train of pulses shown in line 5 of Figure 8. 

The leading edge of the second pulse in this train will be two clock periods from 
the corresponding edge of the first pulse, the third pulse will be only one period 
from the second, the fourth again two periods from the third, and so on. Note 
that the time difference between the last of the five pulses and the first of the next 
train will be three clock periods. During this part of the pulse train, the ENABLE 
gate to succeeding stages goes LOW, and a pulse may be inserted by one of the 
following DRM's. 

D. Cascaded DRM Simulation. The computer was programmed to keep 
track of the number of oscillator pulses which occur between output edges of a 
series of cascaded DRM's. Summing the number of occurrences of each gives 
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Figure 8. Output Waveforms from the SN 74167 DRM. 
Program inputs are shown on the right. 
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Figure 9. 


Observed Outputs of the SN 74167 DRM. Compare with lines 2 , 4 , and 5 
of Figure 3 . The clock input is shown at the top of each photograph. 
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information which can be translated into a mean -squared -error figure for phase 
fitter. The program SYN, which simulates the DRM chain, when programmed to 
multiply by ,52224, is given in Appendix C. 

Since the first DRM will determine worst-case output in terms of phase 
fitter, it is not surprising that in this case there are three possible time differences 
between output pulses. In terms of the oscillator period, the time between outputs 
is one, two or three periods. These results are summarized in Figure TO. 

E. Realization of Omega Frequencies. 

1 , Post-division by 4, In order to obtain a frequency of 652.8 KHz 
from a 2.6112 MHz source, division by 4 is required. This division can easily 
be performed by a pair of flip-flops. Using the SYN program with post-division 
by four (see Appendix C), it was determined that four distinct time differences 
between output pulse edges are produced, as seen in Figure TT . Calculating the 
weighted mean of these pulse occurrences in time gives: 


X = (output pulses) X (number of clock pulses between output edges) 

g output pulses 


= (4448 X 6) + (16672 x 7) -t (23328 x 8) + (7776 x 9) 

52224 


- 7.65931 clock pulses. 


As a check. 


f.^ (clock pulses/second) 


X (clock pulses) 


= f. 


out 


or. 


5 X 10^ 
7.65931 


652.8 KHz (=10.2 KHz X 64) 


as expected . 
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3 




2 

Weighted Mean 1.91483 Clock Pulses 


Figure 10. Output of the DRM Synthesizer with No Post- 
Division Circuitry, In Terms of Input Clock 
Pulses. Baseline numbers are number of input 
clock pulses between output pulse edges. The 
lines represent frequency of occurrence of 
output pulses. 


- 18 - 





6 


7 ij 8 9 

Weighted Mean 7.6593] Clock Pulses 


Figure 11 . Output of Synthesizer with Post-Division by 4. 
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The weighted mean can now be used to calculate mean-squared error. 


€ 


Y^/tlme difference between pulse 
u (.occurrences and weighted mean; 


X (number of pulses) 


2 output pulses 


=|^[(7.65931 -7.00000)^]^ x 76672 
+ [(7.65931 - 6.00000) r]^ x 4448 
+ [(8.00000 - 7.65931 )rf x 23328 
+ [(9.00000 - 7.65931 )r]^ x 7776 J ^ 52224 
= ,69117 'P- 


where: 


r 


f. 


in 


For a 5 MHz clock. 


_2 

& = .69277 


1 

5x 10*^ 


-14 2 

= 2.77x 10 '^sec.^ 


2. Post-division by 8. Since the mean-squared error figure is a 
function of the clock period, it is a reasonable hypothesis to suppose that reducing 
the period will result in reducing the mean-squared error. To test this hypothesis, 
the SYN program was modified to perform post-division by 8 (see Appendix C), 
and a 10 MHz clock was assumed. The resulting output, illustrated in Figure 12, 
shows five possible positions for the resulting output pulse in terms of the input 
clock period. The weighted mean of the output is 15.31857 clock pulses, which 
yields an average frequency of 

10 MHz 


15.31857 


= 652.8 KHz, 



13 


14 


151^16 17 

Weighted Mean 15.31857 Clock Pulses 


Figure 12. Output of DRM Synthesizer with Post-Division by 8. 
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as expected. Mean-squared error In this case is 3.12 x 10 sec. , an 

improvement of an order of magnitude over the divide-by-four case. 

3. Pre-division by 4 . As a point of reference, it is instructive to look 

at the output produced by the DRM chain when driven by a 1 .25 MHz clock, or 

a 5 MHz TCXO divided by 4 (See Figure 10). The weighted mean in this 

-13 

case is 1 .91483 clock pulses, resulting in a mean-squared error of 2.41 x 10 sec." 
F . Conclusions. In the case of the DRM synthesizer used in. the Omega 

/ 

receiver prototype at Ohio University, this study has shown that significant 
reduction in phase jitter, and a corresponding reduction in the error introduced 
into the system, is possible by doublirig the TCXO frequency at the input and 
dividing by 2 at the output. These results are summarized in Table 1 . 



f- 

[ m 

1 

Rate 

Multiplication 

Jl 

e 

1 

5 MHz 

f 4 X .52224 

-13 2 

2.41 X 10 sec 

1! 

5 MHz 

I 

• X .52224 f 4 

-14 2 

2.77x10 sec^ 

Ill 

10 MHz 

1 

X . 52224 7 8 

3.12x 10"^^sec^ 


Percent 

Improvement 


I- Il 88.5% 

HI! 98.7% 

I I - 111 88.7% 


(Percent Improvement — 



X 100%) . 


Table T. Phase Jitter Produced by Three DRM Synthesizer Configurations. 


Dividing the frequency of the DRM synthesizer output by powers of 2 
results in successively larger numbers of possible positions for the output pulse 
edges in terms of TCXO clock periods, it is difficult, therefore, to prove a 
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general theorem relating circuit configuration to resulting phase jitter. However^, 
for a given system, the following rules of thumb for programming the DRM's can 
be observed: 

1 , Remember that the first stage determines worst-case phase 
jitter. Therefore, pick as high a programming number as 
possible for the first stage. 

2. Successive doubling of the input frequency and corresponding 
division by two at the output of the DRM chain has been 
observed to reduce phase jitter. 

in an application where reduction of phase jitter is critical, simulation of the system 
is not difficult using the techniques described in this paper, 

G. Recommendations . Frequency synthesizers utilizing rate multipliers 
hold great promise for use in digital systems where the inherent phase jitter can be 
tolerated. Applications in digital phase- and frequency-lock loops appear 
to be especially useful . When programmable versions like the SN 74167 are 
employed, added versatility is possible. For example, if it is desired to navigate 
with the 13.6 KHz frequency in the Omega system, a simple programming change 
in the synthesizer will yield the proper clock rate for driving the DPLL. Noting 
that the ratio of 10.2 KHz to 13.6 KHz Is 0.75, one could program the rate 
•multipliers to 52224/0.75 = 69632 and arrive at the desired 64 x 13.6 KHz rate. 
(The 11 .333 KHz frequency, however, is not so easily accessable.) 

Care must be taken when using the synthesized frequencies derived from 
SRM chains in analog applications. For example, since the synthesizer described 

. . 5 . 

here is periodic in 10 Input clock pulses (or 10 msec, for a 10 MHz clock), only 
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harmonics of 100 Hz will be present in the output. Changing the phase fitter will 
alter the relative power contained in these harmonics. 

IV GENERAL CONCLUSIONS 

A. Phase Jitter as Noise, The graphic illustrations of phase jitter shown 
in Figures 10, 11, and 12 can be interpreted as probabilit/ density functions (pdf’s) 
when the height of each line, i.e,, the number of pulse edges occurring at each 
position, is divided by the total number of pulses. The probability that a pulse 
edge will occur within the limits of the discrete envelope then will be unify. Al- 
though the phase process is not truly random at the output of the synthesizer, such 
a density function will describe the position of the output pulses when observed 
at random instants of time. 

In the Omega application described in this paper, this density function for 
phase noise can be referred to input since the phase detector used in the DPLL sees 
only relative phase difference between the received signal and the reference 
oscillator. Theoretically it then would be possible to obtain a signal -to-noise 
figure with phase jitter from the frequency synthesizer treated as added noise on 
the incoming Omega signal. Resulting system degradation could then be translated 
into a lower signal -to-noise ratio (SNR) in the signal off-the-air. This approach 
toward analysis has several shortcomings, as discussed below. 

Since the signal off-the-air is band-pass filtered, we can again apply the 
theory of noise in narrowband systems with added sinusoids to describe the prob- 
ability density function for phase. In particular, when the noise input is white, 

( 12 ) 

it can be shown' ' that: 
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2 2 


Figure ]3. Phase pdf for sinusoid-plus-white noise in a narrowband system. 



s = power signal -to-noise ratio (SNR). 

This density function is sketched in Figure 13 for © and increasing 
values of s. Note that when the signal vanishes, the phase is uniformly dis- 
tributed over the interval 0 to 2'tf . As the noise level approaches zero, q(9) 
approaches the delta function at Q — 'If . 

If one wants to speak of the contamination caused by synthesizer phase 
jitter in terms of a lowered apparent SNR, he must fit the empirically-obtained 
pdf to the pdf of equation (9). Several immediate difficulties preclude this 
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approach . 

Firsf, the empirical phase pdf is discrete, not continuous. In the 
general case, even the contour of its envelope need not match any curve described 
by equation (9). Therefore, since equation (9) was derived for band-limited 
white noise (the usual case), we would assume that a phase pdf of the discrete 
form had been produced by some other noise source. 

It is apparent that any signal -to-noise ratio which would result from 
relating equation (9) to one of the empirical pdf's would have a meaning differ- 
ent from the one familiar to the design engineer. Such a SNR would describe a 
synthetic pseudo- noise source with no counterpart in atmospheric channel noise. 

For this reason it would be more useful to formulate a figure-of-merit for a 
synthesizer with given empirical phase pdf. One would then attempt, through 
careful design, to optimize this figure of merit. 

In digital applications of rate multipliers, it is desired that the 
empirically-determined pdf be distributed over as narrow a phase interval as 
possible. The MSE value of deviation from the mean, when used as a figure of 
merit, has the advantage of penalizing large phase excursions more than small 
ones. When there is no phase jitter, the MSE is zero and the density function 
becomes the delta function at the mean. Adopting the MSE as the gauge with 
which to measure phase jitter, the designer should attempt to minimize MSE. 

B . Phase Jitter Relative to Omega Applications. How great an effect 
the MSE has in degrading the performance of a given system depends on the par- 
ticular application. In the case of the digital phas^ lock loop application for 
• which the synthesizer design was developed, it is clear that phase jitter will have’ 
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no effect on fhe system while the Omega signal phase lies outside the inten/al 
embroced by the pdf of the synthesizer. Once the pdf overlaps the Omega signal 
zero crossings, however, the phase detector will begin to give erroneous out- 
puts to the tracking filter. The effect of this phenomenon will be to increase 
the lockup time of the DPLL, making it desirable to minimize MSE in the syn- 
thesizer. 

After the loop has attained lockup, the phase Jitter will drive the loop 
out of lock if the pdf is not symmetrically centered about the lock point, or mean. 
For example, after lock is attained in the pdf of Figure 12, the position of 
Omega phase will coincide with the weighted mean. If we consider the pdf to 
be the result of a random process, then 68% of the time the clock pulses will arrive 
too soon at the phase detector, and 32% of the time they will arrive too late. 

On the average, then, the loop will be driven out of lock twice as often in one 
direction as in the other. 

To completely determine the parameters of this process would require 
a dynamic mating of the SYN program and the CORDET simulation of. Part 11. 
However, the present study reveals that the error will be no more than the least 
significant bit in the bi-directional count register. This is true since the interval 
spanned by the synthesizer phase pdf is less than 1 .53 ^sec, which is the minimum 
time interval the phase can be advanced or retarded at the comparator output. 

The net effect of the phase Jitter in the present system therefore. Is to introduce 
a possible one-bit uncertainty. 

This least significant bit represents 1/64 of a 10.2 KHz lane. Examination 
of Omega lane geometry in North America reveals that in no case is it necessary 
to use station pairs which produce lanes of more than 30 nm separation 
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error of 1 /64 lane is the worst-case equivalent of less than 0.5 nm. When the 
major emphasis is on producing a low-cost receiver, this error can be considered 
acceptable until such time as TCXO's are commonly available in frequencies 
compatabie with Omega. 
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A. COMPUTER PROGRAM CORDET. 


1 . General Description. 

CORDET is a logic simulation of the all-digital phase-lock loop (DPLL) 
presently employed in the Omega receiver prototype-. 

Output provides graphic display of the timing diagrams of the circuit, 
plots of phase difference vs. time, and a histogram of phase difference between 
the DPLL output and the object-lock frequency. 

Provisions have been made to vary the number of bits of integration, the 
signal -to-noise ratio of the input signal, and the initial phase difference between 
the DPLL output and the object-lock frequency. 

2. Capabilities and Limi tations. 

a. The program assumes the presence of a local oscillator in the circuit 
with a frequency of 64 times the object-lock frequency. 

b. • Compile and load time is 15 seconds. Execute time is approximately 
0.15 seconds per cycle of object-lock frequency specified (when only phase differ- 
ence and histogram are plotted). 

3 

c. Approximately 21 x 10 bytes of storage are used. 

d. Subroutines HISTG and RANDNR are called from the Ohio University 
computer center request-call library. RANDNR uses the multiplicative congruential 
method of pseudo-random number generation. HISTG is a plotting subroutine. 

3. - Technical Description. 

Language; FORTRAN iV 

Number of arguments: 1 1 

1. Data card 1: NCYC, MNUMB, NNUMB, 

DCYC, DELPHI, NOISE, A 
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Mode: 


INTEGER (1,2, 3, 6,), 
REAL (4,5,7) 


2. Datacard2: TD, PHA, HiSTOG 

Mode: INTEGER 

3. Data card 3: TITLE 

Mode: ALPHANUMERIC 

Argument definitions: 

NCYC - is the number of object-frequency cycles to be 
considered . 

MNUMB - is the number of bits in the up-counter. 

NNUMB - is the total number of bits in the up-down counter. 

DCYC - is the duty cycle of the m.onostable. 

DELPHI - is the initial phase difference in degrees of lag 
(DPLL output lags input). 

NOISE - is zero for no-noise case, 1 otherwise. 

A - is voltage signal -to-noise ratio in rational number form 

(not db). 

TD - is 1 if timing diagram are desired, £ otherwise. 

PHA - is 1 if phase difference vs. time plots are desired, 

0 otherwise. 

HiSTOG - is! if histogram is desired, 0 otherwise. 

TITLE - is title for histogram plot. 

Figure A-1 represents the positive-logic implementation of the DPLL. 
CORDET performs one pass through a DO-loop for each pulse of the 
64 X 10.2 KHz clock. New values for the logical variables indicated in the dia- 
gram are assigned on each pass. Design changes In the circuit will, of course, 
necessitate changes in the logical statements in CORDET's main program. 
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Figure A-1 . PosiHve-Logic ImplementaHon of the DPLL. 


4. Program Listing and Sample Outputs. 
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INITIATE ANO DIMENSION- 

INTEGER^2 REF, 8( 24 J , INOI TD^PHA , HIS TOG 
INTEGER PHADIF ,PH( 650) 

LOGICAL LMONO,UP, I NDL QC , U P P E * 0 N , L DM EG A , £ N0» LMONO 2 » EX ^ POMEGA 

OM AND AM MUST BE OIMENSIONED TO 2* { 2 * *MNUM B ) 

LOGlCAL^l 0M( 129) ♦ AM{ 129 ) 

DIMENSION ARRAY (650) 

EX = .FALSE, 

LM0N02 = .FALSE. ’ 

INDLOC = .FALSE. 

POMEGA ^ .FALSE. 

READ (1.13) NCYC , M MU MB , NNUMS . OCY C , DELPHI ,NOISE» A 
FORMAT (3I5.2F10 .2 , t 2 ,F5.2 ) " : 

READ (1,14) TD ,PHA ,H ISTOG 

format (315) 

THETA IS INITIAL PHASE SHI FT I «DELPH I /360 

theta = DELPHI >^0.01745329 * ' / ' 

SINE = 5IN(THETA) 

COSINE = COS (THETA) 

00 57 1=1 , 129 
OM< T ) -.FALSE.' ’ * 

AM( I FALSE. . 

CONT INUE 
JJ = 0 
JJ J = 0 

II = 0 
MONO = 0 
M = 0 

N = 0 

MREGM - 2*^MNUMB 
MREGN = 

HALF = FLOAT (MREGM )/2. 

MHALF = HALF 
K = MREGM 

NCYCP = NCYC^MREGM 
NCYCT = NCYCi«=. 1 
X = ( DELPHI/360. ) ^MREGM 
IX = X 

III = 0 
JII = 64 

TNDCH = 0 * ' 

OO-LOOP BEGINS WITH OCCURRENCE OF FIRST CLOCK PULSE- 

DO 12 I=1,NCYCP 

M IS NUMBER STORED IN UP-COUNTER- 

M = M+1 

IF { M-MREGM) 10,11,11 

M - 0 

CONTINUE 

BINARY NUMBER IN UP-DN COUNTER IS TRUNCATED DUE TD SIT 

INTEGRATION. NTRUN IS RETURNED 

CALL TRUN ( NNUMB , MNUMB , N , N TR UN , B ) 

IF M=NTRUN, REF GOES HIGH 

IF (M-NTRUN) 1,2,1 

REF = 0 

GO TO 3 

REF = 1 

CONTINUE 

MONOS PROVIDES OUTPUT OF MONOSTABLE. K IS TIMER FOR MONO 

LMONO IS LOGICAL OUTPUT OF MONOSTABLE. 

K = K+1 

CALL MONOS ( K , R FF , DCYC , MRE GM , L MONO > . 

LOMEGA IS OMEGA SIGNAL ZERO CROSSINGS. 

A IS VOLTAGE SNR. POMEGA IS TRUE PHASE OF OMEGA SIGNAL. 

IF ( IX.EQ.M) GO TC 30 
POMEGA - .FALSE. 


GO TO 31 
POMEGA = .TRUE. 


ORIGINAL PAGE IS 
OF POOR QUALM 


CONf INIJF 

IF (NOISE-EO.O) GO TO 32 
IF (POMEGA) GO TO 34 
GO TO 51 

INDCH - -14INOCH+I 
TF ( INDCH, EQ. 1 > GO TO 37 
00 33 IJ'< = 1 , 129 

AM(IJK) =. FALSE. 

CONTINUE 

NOISE IS ADDED TO QUADRATURE COMPONENTS OF. O, MEGA SIGNAL. - 

CALL ANORM (AN) 

ANUM = A^f'SINE + AN 
CALL ANORM (AN) 
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OENOM = A*CGSINE+AN 
THETAT = ATAN2 {ANUM,DENOM^ 

PCYC = ( THETAI/6.2831 8) ♦N'PEGM 
INDEX = PCYC 

IF ( INDCH.EQ, 1 r GO ‘TO 50" 

AM- ( MREGM+INDEX) =.TRUE. 

GO TO 51 

50 OM (MREGM+INDEX) =.TSU£« 

GO TO 51 

37 DO 52 IJK=1 ,129 

OM(IJK) =. FALSE. 

52 CONTINUE 

GO TO 39 ^ 

51 CONTINUE 

JI I = J I I+l 
III = III’+l 

IF (III. EQ. 129) III - 1. 

IF( JII.EQ. 129) JII =1 
LOMEGA = OM{ I I n ^OR. AM ( Jl I ) 

IF { { I .EQ.M) .AND.{ POMEGA) ) LOMEGA = POMEGA 
32 IF ( NOISE. EQ .O’) LOMEGA = POMEGA 

C LMONO AND LOMEGA ENTER PHASE DETECTOR- 

UP = LMONO. AND. LOMEGA 
DN ={.NOT. LMONO) .AND. LOMEGA 

C PHASE' DETECTOR OUTPUT FEEDS UP-DN COUNTER " " 

IF (UP) - N = N+1 
IF-(ON) N = N-1 

C COUNT REGISTER MUST CONTAIN POSITIVE NUMBER 

IF(ND 8,9,9’' ■ *' " 

8 N = MREGN+N 

9 CONTINUE 

C UP-DN. COUNTER COUNTS MOOULO-MREGN- 

IF IMREGN-N) 23,24,24 ■ .... 

23 N=0 • 

24 CONTINUE 

IF (TD.EQ.O) GO TO 501 
IF ( I-l ) AO, 40, 70 ■ 

40 CONTINUE 

C TIMING DIAGRAM IS PLOTTED— 

PRINT 500 

50 0 FORM ATI ’ 1 » , 2X , » I • , 3X , »N* ,3X , » NTRUN* ,3X,*M*', 4X,»REF» ,3X;» LOMEGA* ,2X' 
, » POMEGA* ,4X , ’ UP ' , 6X , * DN » , 4X , ’ LMONO » , 3X, »LM0N02* » 2X, • INDLOC * ,4X, *E 
^=X‘//) 

70 CONTINUE 

CALL PLOTl { T ,N,NTRUN,M, REF, LOMFGA , LMONO ,LM0N02 , ON , UP , INDLOC'»EX,PQ 
ft.MFGA ) 

501 CONTINUE 

IF { (PHA .EQ.O ) .AND . ( HISTOG.EQ.O ) > GO TO 12 
IF'{ .NOT. POMEGA) GO TO 65 
JJ = JJ+1 

IF ( JJ .EQ. 1 ) GO TO 28 
AJ = JJ/10. 

IJ = AJ ■ 

IF (IJ.N5.AJ) GO TO 65 
28 CONTINUE 

IF (M.LT.MHALF) GO TO 20 
IF { NTRUN. LT..M HALF) GO TO- 21 
22 PHAOIF = 2*(NTRUN-M) , 

GO TO 45 

21 PHADIF = (NTRUN + MREGM-M) *2 

GO TO 45 

20 IF { NTRUN, LT.MHALF) GO TO 22 

PHADIF = 2^( NTRUN-MPEGM-M ) 

45 CONTINUE 

C ARRAY IS LOADED FOP PLOT SU3R0UT INE- 

JJJ = JJJ+1 
PH(JJJ) = PHADIF 
65 CONTINUE 

12 CONTINUE 

IF (PFA.EO.O) GO TO 502 

CALL PLOT (PH.JJJ) 

502 CONTINUE 
IF (HISTOG.EQ.O) GO TO 503 

A.RRAY IS LOADED FOR HISTG SUBROUTINE 

DO 27 I = l.JJJ 
ARRAY (I) = FLOAT (PHI I)) 

27 CONTINUE 

DIMENSION TITLE(20) 

READ (1,146) title 
FORMAT (20A4) 


^ Poor ^xjauty, 


146 
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* I 1 
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10 

‘ 11 * 
12 

19 * 


14 
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1 5 

1 6 
1 7 


DATA AST/^^^/ 

histogram xs plotted 

CALL HIST6 < J J J , A Y t T I TLE ♦ i » AST } 

CONTINUE ^ ^ . 

STOP 

SUBROUTINE TRUH < NNUMS » MNU M3 » N , NT RUN f B } 

THIS SUBROUTINE CONVERTS DECIMAL INTEGERS TO BINARY FORM 
AND TRUNCATESi. SHIFTS DOV/N ‘REMA IN ING SITS AND RETURNS " 
INTEGER^^^a 8<N»NUN8) 

00 9 T ^ 1,NNUMB 
I) ^0 

CONTINUE 

I=sO 

DUM = N 
OUM-DUM/2^ 

r = I*fl ' 

IDUM==DUM 

IF (DUN-IDUM) 2t3,2 
8 ( 1 ) == 1 
DU?4=:DUM- 
IF (DUM) 4,4,1 
8 ( 1 ) = 0 
IF (DUM) 4,4,1 
CONTINUE' ■ ' ^ 

THE BINARY NUMBER IS TRUNCATED 

DO 11 r ^ 1,MNUMB 

B(I) = 3( !+( NNUNS-MNUM8) ) 

CONTINUE ■ ' ■ 

THE TRUNCATED BINARY NUMBER IS CONVERTED BACK TO DECIMAL 
NT RUN = O , ■ ' 

DO 10 T^l^MNUMS 

K^BCI )^( 2-^( I--1) 1 

NT RUN ^ NT RUN -F K 

CONTINUE 

RETURN 

END 


SUBROUTINE MONQS ( K , , CCYCir MREGM , LMONO 1 

LOGICAL LMONO 
INTHGER^2 REF 
X - OCYC^MRHGM - 1 * 

I X “ X 

IF iREF) 10^10, T2 
IF (K-IX) 14*14,11 
MONO = 0 * ‘ 

GO TO IS 

IF (K--IX) 14,14*19 

MONO - 1 

K-0 

GO TO 13 
MONO 1 

IP (MONO) 15 *15 ,1 6 
LMONO ^ *FALSE* 

GO TO 17 
LMONO = ,TRUE. 

RETURN 

END 


OP POOR QUAim. 
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SUBROUTINE PLOT tPH,JJJ) 

PEAL L INEI 129 ) 

INTEGER PHADIF,PH{ JJJ ). 

DATA DOT/* «•/, PLUS/* +*/tBLK/» •/»VLlNE/‘|’/*AST/**»/- 

DO 12 JJ = 1 , JJJ 
IF ( JJ ,NEc 1 ) GO TO 5 
PRINT 8 

8 FORMAT { M * »45X p ’PHASE DIFFERENCE TM DEGREES VS. TIME IN MSEC.* 3 

4 PRINT 2 

2 FORMAT (3X. '-180 ' ,28X» «-g0' .30X. » 0 ’ »31 X t *90* ,28X. • 180* ) 

6 PRINT 7 

7 FORMAT (4X»« + ‘,16<* +*)) ‘ 

5. CONTINUE 

NSP = PH(JJ) + 65 
IF <NSP.GT.129) GO TO 15 
IF ( NSP.LT.l >• GO TO IS ’ 

GO TO 3 

15 CONTINUE 

PRINT 23,PH{JJ) 

23 FORMAT (100X«*PH= 110) 

RETURN 

3 CONTINUE 

C CLEAR LINE 


DO 1 J = 1,129 

LINE(J) = BLK 

I CONT INUE 
LINEIl ) ■= VLINE 
LINE(129) = VLINE 
AI = JJ/10. 

II = AI 

IF(AI.NE.II) GO TO 9 
L1NEU7) = DOT 
LINE(33) = PLUS 
LINEI49) = DOT 
LINE <65) = PLUS 
LINE <Sn = DOT 
LINE <97) = PLUS 

LINE <113) = DOT 
9 CONTINUE 

To - LINE (NSP ) = AST 

IF ( AI .EQ.I I ) GO TO 20 
WR I TE ( 3 , 1 1 ) LINE ' 

II FORMAT (4X,129A1) 

GO TO 12 

20 TIME =(FL0AT( JJ) *. C098039216) 
WRITE (3,21) LINE 

21 FORMAT (4X,129A1) 

WRITE (3,22) TIME 

22- format ( *+• ,4X, E8.2) 

12 CONTINUE 

RETURN 
FNO 
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SUa=JOUTINE PLOTl ( I , N,NTRUN»M »R6F.UCMEGA. LMONO ♦LM0N02 , ON » UP , INOLOC , 
^eX .POMEGA) 

INTEGER^a REF 

LOG I CAL ' LMONO ,UP , INOLOC «DN, LOMEGA * LMON02 , EX , POME G A 
INTEGER LINE(66) 

INTEGER 0LK/’ ‘ / , VL I NE/ * i * 

DO 1 J = I ,66 

• LINE! J ) = 3LK ' ' 

•1 CONTINUE 

LINEU > = VLINE 
.00 2 K - I ,S 

LINE (S^K) = VLINE ' ' 

2 CONTINUE 

IF tREF.NE.l } GO TO 3 
LINE(3> = VLINE 

• LINE m = BLK 

3 , IF { .NOTeLOMEGA) GO TO 4 

LINE no) = VLINE 
LINE (8) = BLK • 

4 IF ( .not cPOMEGA) GO' TO 5 

LINE(18)» VLINE 

LINE (16)= BLK 

5 IFI.NOT.UP) GO T0‘ 6 

' • LINE<26)= VLINE . ^ • 

LINE (24)= BLK 

6 r‘F( -NOT. DN) GO TO 7 
L1NE{34}= VLINE 
LINE (32)= 3LK 

7 IF (-NOT-L.MONO GO TO 8 

LI NEC 42)= VLINE 

LINE (40)= BLK 

■ 8 — IF ( .N0T.Liv)aN02) GO TO 9 

LINE(50)= VLINE 
LINE {4S)=.8LK 

9 IF C .NOT. INDLOC) GO TC 10 
L1NE(58)= VLIMH 

LINE (56)= BLK 

10 . IF C.NOT.FX) .GO TO 11 

LIN£{66)= VLINE 
LINE (64)= 0LK 

11 WR1TEC3,12) I ,N, NTRUN.M.LI NE 

12 FORMAT { 1X,2I4 ,16, I5,SX-,66A1 ) 

RETURN 

END 

100 6 8 0.50 0‘. O 1 3,60 

Oil 

PHASE DIFFERENCE IN 64THS OF A CYCLE 


OEiaiNAL PAOB B 
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a. DRM SIMULATION- 
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VI-PSITY □ U I J I SYSTEM 


16: 17:56 


DFC=iM3EP 9, 1« 


;«AMMhR 
)VIC P 


ACCwi^NT 
NA«vi 3 


iOENT 


60 


PcrcCPDS 

SCO 


CQMP 

F 


LIST 

Y 


zX^C 

F 


J^-r-“rr — ^<*’3rPL. ^ _ 

FOUl V4LCNCE (CK,TA tT«VTC %T01 , <f^ENAei-Xl *KU I 1 t J1 »GI T* 

CL*F3 ♦? t »CC »CO) , ( >T''>,r^*tE^»PPC,PPDJ ♦ (NCKt H2 ;? » D2 ♦ A2 ) »tPID,on » 
iC »C 1 ) * ( Fit* • »l ) f ( A *A1 )»( OA ,N1 ),( C«f*N2> , (35^N3) » (A5,N4) ♦ 

Z,Ol ) * (UN! CA^02) » t.G7» Y ) * ( l 2 , 03 , OCa A» ) »( 0A,K2, 1 2^^031 


«CQB* i3^E3) ♦ C GC ♦CS * ♦ A? » J2 ) ^ ( T A , ) , < J3 , H2 > ♦ ( H3* GO , iC^^COU 

^<OD»?/3) , {OCBAP ,.«A),{F3,CA),te3,C3)»( H4 tOABAP ) , ( V4,6NCUT ) 

A * ( <3 ,CB) 1. (L3,Gft ) 

„ LUG 1 C A L =*^ 1 , PW A ♦ OQ 0 ♦ OOr AO , K CK f FM A 8L £. , CK , T A 3 l JC t XP I x.^UJL 


,GI ♦CL*Fl*02*£l »CC*CD*ST9>F2*S2^PPCfP.^D,NCK,C2»82*A2,«lD*Dl^ 
’^cPICiCl = l Al •DA*N'i tCA,N2t55*N‘3*AStNA,N5,2,Cl ,UNICA,02t03* 

9 Lk *D3. CC6A^-i*QA ,02. I 2 . C 3 t OH i 13 . 03 .QC* G2 . ‘^2 ♦ A3 ♦ 32 ♦ T4 ^Hl , J3, H2 1 H3 * 
.JGOG3,GC,0r>^,y3,00nAP, Aa,F3»CA*E3»C3TR4,OABA« ,ENCUT tK^tLa s K3,GA_ 

INxrCEP LINE (66) 

U^JICA = •TRUE* 

QA5 *FALSF2 


03 = •FALSE. 

GC =: .false. 

Qo ~ .false. 

QA3AP ^ , . TPUg. 
CCEAO ==: .TPUO. 

ODBAP = .TRUE. 
P4A « .FALSE. 
PRo .FALSE. 


REAP (Itl) EMAfLe.STO.CL 

1 FOR^n^AT (7C4XL1)) 

-C—3”:: DCT^oMiNif INPUTS , 

NENAOL = -KOT.eNAaLF 
DU 6 U=1 *10 

READ (1.7) R ID.RIC.RIB.RI A 

format -C A (AXLl >) - 

on u 1 I ,AC ' " ‘ 

F3 =: F1.CR.F2 
E3 - Fl.CR.F2 

L3 ^ L1.ANO.L1L _ ‘ _______ 

K3 <l.AVn.<2 
T4 :=i I 1 . AND. I 2. AND. 13 
J i ^ 3 1 .and. 32 
H3 ^ HI ,CR.H2 

— ™GENtRAT.E CLCCk'P'JLSF " 

ICK - --l^lCK + l 

IF ( ICK.FG.l ) GC TO 2 

CK_^ .FALSE. 

GO TO 3 

2 CK =: .TPUF. 

3 CONTINUE 

,Crj:3rrc.*^-rLEA0 .Epgf .0F_CK^ to toGEPS., F FLO PS - 

CALL TFF (CA.PpA.GA. tA .QA, CAPAR) 

CALL TF^ (CD>PRF.G==»iTO.OS.Ct?BAR) 

CALL TF- CCC.PFC*GC.TC.OC.GCaAP) 

_CALL,^ TFF (CO.Pf^C ,GC>TD.j:5t: ,CDi^AJRJ 


OetFOMtME: CUTOUT 

NCK ^ .\0T.CK 

^4 o .fjCT . (Mi * AND. iVo. AHO.?43> 

AND .JD,3 


C4 ^ C I . Aun.CS .AND.C3 
85 =:■ 31 . AND.«i>.AND.03.ANC.B4 
AS ^ Al . ANO, 42 . AND .A3 . AMC. A4 
NS a ,, . NO T . ( .N I . 3R . N 2 . CR . N 3 . CR . N4 


03 = ( .tiCT.COi .AN0.02) > 

IP (l.FQ.l) PP|NT7 .Pjr.PICfRIBfPI A 
IF < I .EC .) / PP TNT 5 

RM AT { 2 3 X . • C K » , 7 C . » V * , ax , t 0 A e., 6 X . «il3 \ . 6X,#/ QQJjt 6 X t.r.0P-lj 

♦EN* . 3X» *FNPUT < ,4X , * 003AP* ) 

CALL PLPTl <I*P ♦ ICK. ICK.CK. V.OA *OB.QC.OD.SNAeLE»ENCUT.OOBAP> 
CONTI NUT 

, _CO^.T tNUS 

STOP 
END 


F.CU1 valence data map 


LOCATION 

OGOORC 

OGOO^D 


SVMriOU 

TA 

LI 


LCCATinM 

t'vOr'tC 


SVMDCL 

TO 

K1 


LCCATIOPi 

CGooec 

COOf'fiO 


SYMBOL 
TC 
I 1 


LOCATION 

OOODOC 

COOCED 


OKIfflNALmEB 

01 K>0R <SDA1®^ 
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DECEMBER < 


IVEPSITV 0 {} 1 j 1 system Fortran iv g-levEl compiler 

SUBP^OOTINE PLOT I < 2 I * J J Tkk »MM^Ll tLRTl 3»LA ,L5 f 

logical* I LI *L2»L5.LA^L5,Lo*L7,LatLP 

INTFGFP trNr{66) 

IMTf;CEP PLK/» VLINE/* J •/ 

00 1 j = 1 ,tva " — — 

LlNLIJ) « PLK 

1 CONT iNUF 

LL^bJ'LI.l )._«L.VL.ItiE 

DO 2 < =s 1^8 
LINE (S^O = VLINC 

2 CONTINUE 

t •NOT.L! > OO TO 3 

LiNfcO) ^'VLINE 

LI NE ( 1 ) -=5 BLK 

3 IF ( •N07 .L2) go to 4 

LINE VLtNg., ' 

LINE <a) 3= BLK » ^ 

4 IF i .NOT»L3) GO TO 5 

LlMEdST* VLINF 

: L I NF ( 1 6 )^^L , 

5 IF (*N0t*L4>'G0 to' 6 ^ 

LINE<26)3!:: VLINC 

* LINE 124)- PLK 

NGJ.r L« L*jGQ_tO_7_: j 

LINE(34)ss VLINE 
LINE I3P. >s: BLK 
7 IF (♦NOT*L6) GO TO a 

L I <42)=? VLINE 

L I NC“ { 40 ) -"IfJLK' — 

6 IF ( ^NOT .L7) GO TO 9 

LINE(50)*x VLINE 

LINS, ?4SJ,F.._0LK ^ 

9 IP (•MCT*LB> GC TO’io* 

LINb(5e>-=5 VLINE 
LINE (56)as 9LK 

10 IF (cNPiT^LQ) GC TO 12 

L tNfc VLINE" ” 

LINE (64)« BLK 

11 WRITE <3*12) LI * JJ *KK*MM .LINE 

12_ format, UA*2I;^ IEj I5^5X,6FAl > 

RETURN. Vy " --- 

END 

SCALAR MAP 


LOCATION 

OOOGE4 

ocooFa 

JOOCIOO __ 
OOOiOB 


SYMBOL LOCATION SYMBOL LOCAT^ION 

VLINE OOGOES J OOOOcC 

KK CrOOFC MM OOCICO 

. L4 ^0 n rj 07 ^L,5 _ OC 0 X 0 i? 

L9 00 0 IOC ' ^ 

ARRAY MAP 


SYMDOL 
< _ 
Li 


LOCATION 
OOOCFO 
0001 C4 


JLSCAXLQt^L 

DOOi 10 


000218 


-SYMS OI L OCATI phL 


.iyM8.oL 


SUBPROGRAMS CALLED 


, LgCA T l.QN, 


,.^YMSOt 


LABEL *MAP 


^.^iJ3S:^jrSQIiL^ 


_JLa.QAJ LQI^ UA£LLL L PJ" aJ: 

0003BC 2 C»C05F2 

0OO4S0 7 0CC49r 

000510 12 00C5SC 

MEMORY REQUIRgM F NTS 000 5S4 BYTES 




. lo^catic 

C0C436 

CCC4BC 


C00444 

0004DA 


ORIGINAL PAGE IS 
OP POOR QUALirYi 
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OHIO UMivrRSiTY O U I J J 


FO??T«?AM IV G-LSVFL CC^^PIL£R 


or 01 
00 0 2 
CC03 
0004 

cro^ 

00 Do 

coa? 
ocos^ 

0009 
CC 10 
.cm 
.P0X2 

CC 13, 

00 S 4 

CCI5 ' 

,00 1 6 

0C17 
CO 18 


SY ^BCL 

PP 

T 


SU«'-^Ol^T INr TFF ( C*PP,G*T« C*CBAR) 
LOGICAU^l C * FP , r »GBA!? , G 
IF ( •NOT •(PQ^Or^.C) ) GC TC 1 

CPP) GC_TC 

O'* »F aL SF « 

GO TC 3 
COKTIK‘U^ 



OaAPst.HOT^O 

IF ( -MOT. (PP .ANn^C ) ) R£TURN 
0 = .TP^*£. 

CSAP=5*TRJE# _ 

RCTUPM ^ ' “■ 

IF i .MOT.G) RETURN 
IF (T) G=.N0T.G 

QRArv-=.MOT*0 _ 

RCTUmN 

FMD 

SCALAR MAP 


LOCAriOM 

OOOCBC 

ooooes 


SYMSOL 

c 


LOCATION 

OOOCBl 


SYMBOL 

0 




label location label location l 

e 0OO1B6 ^ 3 COOICC 

TOTAL MEMORY RFMfNTS 000364 BYTES TFF 

/DATA 

module acdr length Entry points <E) 

MAlN44i^ OFOOO 0C6A4 E-MAIN^^. OFOOG X-Il 

.-ELctj aEaA.e. .C050_4 e-plchl ofsa^ 

TFF## 0FC6C 002^4 E-TFP 0FC6C 

J^COLLE successfully LCACED 300DO BYTES QF ST! 


LOCATION 

00,0083 


LOCATION 

000208 


SYM80L 

OSAR 


Entry points <E). external cx) ano common refe 

OFOOG X-ISCOM^^ X-TFF 

_ 0F6AS 

0FC6C 

30ODO BYTES QF STORAGE REMAINING / EXECUTION BEGUN 



ORIGINAI. PAGE IS 
Ot POOE QUALIIY 
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C. COMPUTER PROGRAM SYN 
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io:is:27 


FeSRUARY lA# 1 


CVERSXTY , Q U I J I SYSTEM 


iRAMMER 
BERLIN K 


ACCOUNT 

N25A6 


XDENT 


TIME 

ICO 


RECORDS 

1500 


COMP 

F 


c- 

c- 


c — 


c- 

4 


c- 

5 


7 

c- 




.8 

10 

2 

9 


LOCATION 

0000D4 

OOOOEa 

OJOCFC 

OOOlOl 


LOCATION 

000JD4 


HISTG 

logical*! SYN1*SYN2# SYN3,SYN4,SYN3» SVN*HNAaLE 
INTEGER PREV 

DIMENSION TITLE ( 20 ) ♦ D I F ( 6000 > » I DIF( 80C0> 

PREV = 0 
Ml =: O 

K “ O 
KK =s 0 

DO 9' I I 1, !0 
N =; 0 

IF— PREV PREV - lOOOO " ” 

'DO 8'ICK ^ !• 10000 

FIRST DRM IS PRCGRAMMEC TO S 

SUCCEEDING STAGES ARE ENABLED WHEN ENABLE ^ *FALSE«*-» 

I = MOO IICK.IO) 

SYNl = I U *EG^2 ) .ORo 1 1 oEG.4> ^OR. ( I .EQ*S) oORc{ I .E0^7) •OR*t I. 
ENABLE r •NEcO) 

SECOND DRM IS PROGRAMMED TO 2- 

IF CENABLE! GO TC 4 • 

M = M+1 
I = MOOCH, 10) 

SYN2 = C ( I ,EG*3} •OP* t I •EC,6) I 
ENABLE I #NE*0) 

THIRD DRM IS PROGRAMMED TO 2— 

IF { ENABLE) GO TO S 
4 = J + 1 

I ^ MODI J,10> ^ “ ” 

SYN 3 = <U*EQ-3)*CR*II*EQ*8I> 

ENABLE =sU*N£*0> 

FOURTH ORH Is PROGRAMMED TO 2 - 

IF I ENABLE) GO TO 6 ^ 

K ^ K+1 

I = MOD (K,IO) 

SYN 4 = {(I*EQ.3)#0R,( l.EQ,a>) 

ENABLE -NE.O) 

FIFTH DRH IS PROGRAMMED TO 4~ 

IF C ENABLE) GO TC T 
KK - KK 4^1 

I = MOO(KK,lO) ’ ' ' 

SYN5 = ( U«.EG.2).CR*CI*£C*4)*OR*(I,EQ#7)*OR.<r*HQ*90J 
CONTINUE 

-CLOCK OUTPUT IS SYNTHESIZED V^ITH OR-GATE - 

SYN = SYM •0R*SYN2*0R*SYN3.0R*SYN4«0R*SYN5 '*■ — 

-^--NQRMALrZED PHASE CIFFERENCE IS CALCULATED 

IF {.N0T*3YN) go TC 3 
N == N^' I 

lOrF(N) = ICK-PR6V 
PREV =: tCK ^ 

SYNl = .FALSE, 

SYN2 .FALSE* 

SYN3'- .FALSE. ' 

SYN4 ^ .FALSE. 

SYN5 = .FALSE* 

CONTINUE 

OO 10*' I =s 1,N 

DIF<I) = FLQATdCtFl ni 

CONTINUE 

DATA AST y»**/ 

READ <i,2) TITLE 
FORMAT C20A4) 

CALL HiSTG CN,0IF,TITLE,1,AST) 

CONTINUE 

STOP 

END 

SCALAR MAP 


SYMBOL 

M 

N 

enable 

SYN 


symbol 

DIF 


LOCATION SYMBOL 

000008 U 

COOOEC ICK 

OOOOFO SYN2 

C00102 

ARRAY MAP 


LOCATION 

000154 


SYMBOL 

IDIF 


LOCATION 

OOOODC 

OOOOFO 

OOOOFE 


LOCATION 

007E54 


SYHEOL 

K 

I 

SYN3 


SYMBOL 


LIST EXEC O 
Y F 


EQ-9) I 


LOCATION 

000060 

C000F4 

OOOOFF 


LOCATION 



normalized frequency spectrum 


PLOT 1 


3182 

♦ 





3139 


* 

■* * - 



3C96 

* 





3053 

* 





3010 

♦ 

+ 

+ + 


+ 

2967 

" * 



" 

■ " 

2924 

♦ 





2881 






283 8 

♦ 





“2795 


4*' 

“ + 



Z75Z 

♦ 





2709 






2666 






2623 

" ‘ ♦ 


.. - 



2580 

* 

+ 

•+ ^ 

. + 


2537 

* 





2494 

* 





2451 




■" * 

— 

2408 

* 





2365 

* 


+ <• 

4- 

4* 

2322 






2279 



" *" — ' “ 

' * 

■ — * — - 

2236 ^ 

♦ 


* 



2193 






2150 

* 

-♦* 

+ + 


. 4 - 

2107 , 

* 





2064 , 






2021 i 

* 





197 3 ; 






1935 ^ 

* 

+ 

' + + 

■t' 

4 * 

1892 

* 





1649 






180 6 

* 


' 



'*1763 



*'■■ * * 


- 

172C 



+ 

4 - 

+ 

1677 

* 





1634 

♦ 





159 1 


" * 





1548 

* 





' 1505 

+ 

4* 

+ + 

+ 

4- 

1462 

* 





1419 



” 



1376 

1 * 





1333 j 

* 





1290 

* 

4- 

+ + 


4- 

1247 






1204 






1161 

** 





1118 

4 + 


, 

'* 


1075 , 


4- 

- 4 . 4 . 


4- 

1032 

*♦ 





989 

♦ 4 





946 






903 




- 

- ’ - 

860 


+ 

'f' 

+ 

4* 

81 7 

** 





774 






731 


* 

— . 

- ’ 

•- 

68 8 * 






645 


+ 

-I- + 

+ 

4* 

602 






559 



" 

. -* 


516 : 






473 

!}E»^ 





430 



4 + 

+ 

4- 

‘ 38 7 






344 

* + ic 





30 1 






258 

**♦ 





21 5 

*** 

+ 

+ 4 - ' 

+ 

4- 

172 

4! + * 





129 

’If** 





86 

*** 
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*♦* 

1 ^ 

1 t 




s 

1 

r 

6 

11 16 

^ 1 

21 

26 


NUMBER 

OF POINTS. 

REPRESENTED 


+ 


■i* 


♦ 


+ 


-f- 


* 


+ 


'I* 


+ 


-F 


- 6 - 






+ 


-I 

31 

IS 


. . , + + 


+ ‘ + “'F 


i* 



*f* + + ■J* ' 


+ ^ + _ + + + 



t ♦ + 





+ 


'j- 


■i* -ft* 


+ 


+ + ■► . + . 


+ + • 4 , 4. 


j I I I |, 

36 41 46 51 56 

5223 


*• 4 * 4 . 


+ “ 4- ~ "V 


.4* + ♦ 




+ + + 


i- + + 


+ + 4. 


4- + + 


+ _ . _ 4* + 


+ + +' 


+ 





+ • + .+ 




4* ^ + 


€»i 66 ?1 . 

NUMBER OF Pf 


ACTUAL PLOTTED DATA LIMITS ARE FROM 1 TO 


INTERVAL FOF 




v=rr?s'iTY n u I J I 


i i : 55 : 1 a * m a» c u 31 1 15 7S 


A-CrUNT IDCNT TIV5 «£COPDS COMP LIST EXEC 0 

K L^fLO 5.9^ ^ ^ 


c 


To 


c 

c 


c 


c 


c- 


-C 

6 


7 


C 


c 


u 
^ ’ 


9 


LOCATION 

300054 

OOOOcB 

000 JEC . . 

0 0 0 1 I 0 


SYNTH'^SIZE^ MOO*=L ^1TH POST-DIVISION BY A 

L DC, I CAL SY^41 ^ SYN2* SYN3 . SYN4 » SYN5 • SYN , CNABLE « OUT 

IVTb'C'^P POPV 

T'JTI lA . 

Dlf-'F-NS ION ICNT <f 0) 

DO 10 I = i.eo 

I_C NTJ I ) = 0 

COTnI INU5 
IDIP = I 
P*^CV = 0 

V = 0 

J = 0 
K = 0 . 

KK = 0 

= t .AciTo"^ ’ ~~ 

first ORM is PPOGRAMMEO to 5 

SUCCEEDING stages ARE ENABLED WHEN ENABLE = #FALS£* ' 

I =: ypO ( ICK 0 ) 

SYNl = ( < I *rO ♦ 2 ) .OR . ( I .EQ .4) .OR.JI*EQoS)*OR#(I*EQ«7).OR.(I.EO#9)) 

ENAbLL =(I.NE.O) 

SECOND OP-M IS -DROGRAf^MEO TO 2 

IF (FNJk«LF) GO T O 4 



I = C M , 1 o > 

SYN£ = ( ( I .EC.5 ) .OR. C I .EG.3) ) 

enabl e =(I>NE ._0J 

Third or.m is programmed to 2 

I- <C.vAr5LC> GO TG S 
J = J + 1 

I - .vOD<J.-lO) 

S^N ? = < (T7E0.3 ) . CR. ( I « E0.3 ) ) 

ENA3L? =(I.NE.3) 

fourth D-J'L is programmed TO 2 

I ^ ( FN ADL^ > GO TO 6 ' 

K = K + 1 

I - YOn (<,I0) 

SYN A = ( ( I . fcO .5) . OR. ( I.EC.S) ) 

ENmHLl =(T.NE.D) 

=^'--^frYT^OR~<^ IVIPROGRAMMEC 'TC 4 

I- (un^;»=«lE) go to 7 

KK := KK +1 

I = M OO ( KK J l_p ) 

3YN5 ~Vi i.CQ.2) .OP.< I. ED. 4). OR. ( I.E0.7).OR.(I.E0.9)) 

CONTINUE 

— CLOCK OUTPUT IS SYNTHESIZED V.ITH OR-GATE 

SYN - 3YMI .OP. 3Y NF.DR.. SY^ jtOg,F^S,YN A. CR . S Y N5 

Q'JTPUY IS DIVIDED 8Y 4 

IF (.NOr.SVN) GO TC II 
NN = NM+ 1 

^»;iN -=i MCDJCNN.a) 

OUT = (MN.EQ\o) 

normalized PHASE DIFFERENCE IS -CALCULATED 

IF (CJUT) IDTF - ICK - FREV 

I ^ ( I CT F ) ^ I CNT__(.LDT ET-^!^ 

IP (ruT) Wjev = ICK 

sYMi - .False. 

SYN2 = .FALSE. ’ * . ' 

SYN5 - . FALSE. 

SYN 4 - .false. 

SYN5 = .false. 

CONTINUL 

C 3NT INUL 

0G“9T“= 

P‘?ir;T 500 . I ♦ ICNT{ I ) 

CO NT INUF 

_ P.JPMAT UXt2 l7 ) : 

'stop 

FNC 

SCALAR MAP 


VrMpOL LOCATION SYmpcl 

lOIF 30000S PREV 

KK POOOEC NN 

SYN^t i'C0!00_.„ SYn3 

Mil ”’ C C 0 1 J 4 OUT 


LOCATION SYMBOL LOCATION 

OOOOOC M OOOOEC 

OOOOFO ICK 0000F4 

000104 SYN4 oooioa__ 

000118 



-51 


O H I C U N I V r ^ s’ I T Y ■ J I J ' p* s Y ST F*? ^ TR AN IV G-L F v'El"^ C OMP l"Cfe R 


symbol 

LOCATI ON 

SYMBOL 

1C NT 

0001 1C 


SYMBOL 

location 

SYMBOL 

IRCOMA^ 

000X54 


LAHEL 

location 

LAHFL 


ARRAY MAP 

C A_T I C«N S X M f O L L OCA T i pN 

SUBPROGRAMS CALLED 

.kQ,C^IjQN SYMROL LaCATlON., 

LABEL MAP 


S_Y>lSO! 


SY-MBOt 


LA35I 

6 

500 


_ LAH£L kAkATION LO.Q^T I ON L AJ4£L A. 1 1 ON LAS 51 

10 ooon ici A- 0C0‘’9S 5 OOQAIO 6 

11 O^IOBOC a OOOSOC 9 00062C * 500 

TOTAL MFmOPy REOuI PEMCnT S 0006BA OVT SS MAIN44 

_-/t?.ATA 

MODULE * ACCR'*' LET-TGf H ENTRv POINTS (5), EXTERNAL <X) ANO” COMMON <CJ Rl 

MAIN44<^ 0=^000 005^4 E-MAIN44 OFOOO X-r0CCM/i^ 

MODULE SUCCESSFULLY LCADED 20900 BYTES CF STORAGE REMAINING / EXECUTION SEGl 




09t4S:05 


LIST 

Y 


Vf^SiTY n u I J T SYST{=‘4 — 0'9t4"S“:05 — ” 

ACCCU^^'»’ IQLNif live RECORDS COMP 

< L2n^o Sjoo 50.0 __ 

C — i:OY Sv:gTHFSlZf*«‘ MODEL VfITH POSt-D I VI SION 3Y 9 

L^SICAL 5YN1 »SYNS, SYN9 *SYN4, SyN5*SYN»ENABLEf OUT 
IN4T6(3r- PWEV 

I ^‘T f‘ ^ P I A ^ 

Oi KS fCN ICKT(SO) 

I “ 1*50 
ICKTil ) - C 



ror»« == 1 

OPSV « 0 
M ^ 0 
^ J ^ 0 
K 0 
KK * ^ 0 
NN - O 

"= 1,^00000 

c ►"'-FIkOT 13 PPCCPAMMPC TO 5“ 

c SUCOnHOTNO STAGES ARE ENABLED WHEN ENAaLE = .FALSE*— 

I - VOD U C<*101 

SYNl » t < I .50.2 J .CR. { I .EC.4) .OH • f I » EQ .5 > .OR . ( I . EQ .7 I .OR • C I . 60^^91 ) 
ENAbLF ^(I.Nf.0) 

C DR?v IS cpcGRAMMEC TO 2— - 

IF (FU4“iLL) GO TC 4 



I » W0P(M,10) 

3YH2 » ( ( I .50. 3) «0P. f I oECp.9) ) 

F!±AMLE =(I^nE.O_) * * 

C fnlr^') ORM IS PRCCRA.VVED TO Z- 

4 IF (EWABLF) GO TO 5 

J ^ J Hh 1 « ‘ 

I a ^^GD{a,lO) 

SYN 3~‘'sTTri .EQ .3 ) I .EG.a ) ) " 

ENABLE «(I.NE.O) 

C FOURTH DRM IS PPCGRAFVEO TO 2™-* 

_S Ir^(PNANL5I CO TO 6 

K - k*Y'^ ' “ 

t MOD cn*10> 

SVM 4 ( < I.EO*3).CR.Il.FC.8U 

_ ENAbtE -ei.NE.O) 

CfTff PRDGRA.VJ^C TO 

IF (SNAPLE) GO TC 7 
KK ^ K< Hhl 

U JiOni Ki^ zJJD i 

SYNS ^ iTl .60»2} .CP. ( ! .EC .4 J .OR* < I .€0.7) .0». { I .EQ.9) ) 

7 CONTINUE 

CLDC>« OUTPUT IS SvMTHHSiZFD WITH OR-GATE 

SYN “= SYNl .CR. SVN? .OR .S^M .OR. SVNa .OR .$YNS 

-:C-~~-OuTpUT IS DIVIDED S Y 

( .NCT.SYN) go TC 11 
NN « NN4- 1 
MM 2x M(JO (NN*e) 

OUT IMN^'eO.O) ^ 

C --NOt-PAL r ZFD RHASE CTFFERFNCF IS C ALCULATEO— ™“~ 

IP (ruT) TOIF TCK - FPSV 

IF (mTl__vCNT C iniF> _?_ ICNT ( TCIF ) 1 ; 

IF (OUT) PREV = ICK 
SYNl 2S .FALSE. 

SYN2 ^ .FALSE. 

S YN3 g .false. ^ " 

SYN4 r: .FALSF“;;^ ^ = : 

SYNS .FALSE. 

II CDNTINUr 

a continue 

OS'o r'- 1*50 
PRINT 5f»n, T * ICNTI I ) 

9 CONTInU** 

^;00 FDRM AT i J X *J> 1 7J 

S TCP”* ' 

END 

SCALAR MAP 


MARCH 31* 197S 
Of 


exEc 

F 


LOCATI orj 
000^04 
0O0*'> *^8 
_0 

OOOl'lO 


SYMACL 

Dir' 

SVfj£ 

l4N 


LOCATION 

''ocooa 

^0 CCEC 

_ocoioo 

Q003 14 
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SYMBOL 
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OH I O' U\''l VPO SI T Y 'O "U T-J - f— rsYSTFlS FOpf'^TN- 


IV G-UevfL' C’OMpfLHR' 


symbol 

1 CKt 


^LUCATI OK 

00 o'l Vc 


SY 


APPAY 

— Jk^A.tXQ.N 


f^AP 

_ SYMBni, 


SU0PPOGR A?/S 


.J-OCATI.ON 

calubd 


Sr^SOL. 


_SYMH OL 
ltiCOyV~ 


_LJCA^i nN 
00 0Tt£4 


S VMU r.L 


^L.CCA T I OM 


^s:£MD,qu 


-UDjQAtim. 


JSYMMt 


LOG ATI ON 
“0 

0005DC 


LAFi^'L 


LABFL MAP 
-LJ-^A T I ON 




MAIN44 


J-aC^JJiDNL 

0004 1 0 

ooo6ac 


OOC'5 9« 

total MrMDF<Y rJ^QUI PE'^TNTS C006B4 YT?S 

^UDULB “ AOn»^ L-NCTH (H>rPxtERNAL”x)”^ 


■ L.ASak. 


MAI N4 A. 


Aon»-> 

OFOOO 


LFNCTH 

006n 


r-MA l^{4 A 


»£NTRy 

oForo 


MODULE SUCCESOPULLy LCADHO 30000 BYTES GF^STORACf PEMAINING 
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